Onset of Antiferromagnetism in heavy fermion metals. 
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There are two views of antiferromagnets. The first proceeds from atomic 
physics, which predicts that atoms with unpaired electrons develop magnetic 
moments. In a solid, the coupling between moments on nearby ions then yields 
antiferromagnetic order at low temperature^. The second, based on the physics 
of electron fluids or 'Fermi liquids', states that Coulomb interactions can drive 
the fluid to adopt a more stable conflguration by developing a spin density 
wave.l'l It is presently unknown which view is appropriate at a 'quantum criti- 
cal point', where the antiferromagnetic transition temperature vanishes&~i. Here 
we describe an atomically local contribution to the magnetic correlations which 
develops in the metal CeCug ^Aua, at the critical gold concentration (xc = 0.1) 
where the magnetic ordering temperature is tuned to zero. This contribution 
implies that a Fermi-liquid destroying spin-localizing transition, unanticipated 
for the spin density wave description, coincides with the antiferromagnetic quan- 
tum critical point 

CeCugi'i is a "heavy fermion" compound, a class of metal formed between actinide or 
rare earth elements and transition or noble metals. 'Heavy' refers to the extremely large 
effective masses of the charge carriers at low temperatures, often hundreds or even thousands 
of times greater than the electron mass. Because of their large effective electron masses and 
proximity to antiferromagnetismEHlil, heavy fermion materials are ideal venues for observing 
the competition between metallic electron band formation and magnetism. The large masses 
imply small bandwidths, and the magnetism is easily tuned by external pressurelll, magnetic 
field and chemical composition. Thus, all of the phenomena associated with the competition 
can be readily accessed at low temperatures where it is much easier to collect clean data 
than for other materials, such as the transition metals or their oxideJii0, which also display 
competition between band formation and magnetism but whose characteristic temperatures 
are much higher. Because of a substantial data base&H and the ability to grow large and 
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highly homogeneous single crystals, we have chosen CeCue-xAu^,. as a test material. The 
quantum critical point (QCP) occurs as gold is substituted for the Cu atoms; when more 
than 0.1 Cu sites per cerium are replacedcJ, the heavy fermion paramagnet gives way to an 
ordered antiferromagnelSill (Fig.|I]a). 

Fig.0 illustrates two competing models of magnetic order in rare earth metals that have 
been widely discussed in the context of heavy fermion systemsiUnl. The fundamental pa- 
rameter is the strength of the hybridization W between the localized f orbitals and the more 
extended s, p and d orbitals. When W is small, the extended orbitals form a metal which is 
weakly coupled to the unfilled, and hence magnetic, f-orbitals. Thus, as T is reduced, the two 
subsystems undergo largely independent evolutions towards Fermi liquid and magnetically 
ordered ground states, a situation found in many elemental rare earths(e.g. gadolinium), 
and their compounds. The volume bounded by the Fermi surface is "small" , containing only 
electrons from the extended orbitals. As we begin to increase W, the magnetic coupling be- 
tween the rare earth ions via the conduction electrons and the magnetic ordering (or "Neel" ) 
temperature, T/v, both increase. On the other hand, if W is large, the material develops a 
Fermi liquid of strongly hybridized carriers with a "large" Fermi volume containing electrons 
from both extended and f orbitals. The kinetic energy of band electrons then overwhelms 
interaction effects, with the result that T/v is suppressed to zero for W beyond a critical 
value W^. 

A key scale for heavy fermion materials is the "Kondo temperature" T^'- the temperature 
below which the magnetic susceptibility saturates and the hybridization of the f-orbitals with 
the other orbitals becomes apparent. In contrast to T/v, the Kondo temperature Tk, increases 
monotonically with W. The central open question concerns how this scale behaves near the 
quantum critical point: in particular, whether it remains finite, (Fig. lb) or whether Tx and 
T/v vanish at the same point (Fig. Ic). The first scenario is required if the local moments are 
quenched at a finite temperature above the quantum critical point: here, local moments do 
not play a role in the physics of the quantum critical point and a spin-density wave model 
must be used.iUliil In the second case, local moments are present at all temperatures down to 
T = at the QCP, so a heavy electron Fermi surface ceases to exist at the quantum critical 
point. The present paper gives comprehensive evidence in favor of the second scenario, 
fig.S(c). 

To discriminate between the local-moment and spin- density- wave visions, we have mea- 
sured the magnetic response function x(q, i?), as a function of temperature (T) and ex- 
ternal magnetic fields (H), using neutron scattering and bulk magnetometry. The func- 
tion x(q, £") describes the Fourier components at frequencies E/h of the magnetization 
(5m(q, t) = x{(l,t)Sho induced by a unit external field pulse 6h{r,t) = 6hQe'^'^'^6(t) with spa- 
tial modulation described by the wavevector q. A particularly simple impulse response is 
an exponential decay, ^(q, t) oc exp[^^t], where Fq/^ is the magnetic decay rate, for which 

x{q,E) = A/{T^~tE), (1) 

where A is a constant. The exponential form is typical for paramagnetic rare-earth insu- 
lators where the magnetization is not conserved. To understand how Fq can depend on 
temperature, recall that the zero-energy response is the static magnetic susceptibility, which 
in paramagnetic insulators follows a Curie- Weiss law x(q? E)\e=o = C/ (T + 9{q)) with Curie 
constant C and q-dependent Weiss temperature 9{q). At q = this is just the uniform sus- 
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ceptibility, which can be directly measured using a magnetometer. Comparison with eq. (HD 
leads to the conclusion that Fq = a(T + 9{q)) where a = A/C. Hence, in paramagnetic 
insulators, the inverse magnetic response function 

X-\ci,E,T) = Xo\E,T) + r(q)/C (2) 

is a sum of a non-local, purely q-dependent piece 6'°(q)/C (here with a = 1) and a local 
(q-independent) driving term Xo'^ — (T — iE / a) / C . Thus T and iE/a are interchangeable 
in expressions for x(q, i?), and at critical wavevectors, i.e. the points in the Brillouin zone 
where ^(q) = 0, x = Xo exhibits E/T scaling, meaning that 

Xo = T-''9iE/T) = E-''jiT/E) (3) 

where j{l/y) = y°'g{y)- In the simple Curie-Weiss example, a = 1 and g{y) = go{y) = 
C/{l-iy) with y = E/aT. 

While ii^/T-scaling at critical wavevectors is well-known and establishedlzllil§§, it does 
not capture the general dependence of x on wavevector, magnetic field, and zero-temperature 
tuning parameter (e.g. pressure or composition) for the quantum phase transition. However, 
inspection of the Curie- Weiss example suggests a set of hitherto untested scaling laws and 
data replotting procedures. The first follows from insertion of the T ^ limit of eq.(|^) into 
eq. d^) . The outcome is that x~^{^y E,T = 0) displays E/9{q) scaling, i.e. 

^(T ^ 0) = e{q)--G{E/e{ci)) (4) 

In the simple Curie- Weiss case G{6) = g{6) (with 6 = E/a9), which is easily understood 
if we regard excursions in wavenumber or quantum control parameter as excursions in an 
effective Curie temperature. Letting E ^ instead of T ^ in eq.(^) and again inserting 
the result into eq.(0) gives a second verifiable relation, 

X-^(q,T) = (T" + %r)/C, (5) 

Eq.(^ implies that the traditional Curie- Weiss plot, where 1/x is a straight line as a function 
of T, can be replaced across a series of experiments where the Weiss temperature is varied 
(e.g. by varying q or composition x) by a plot where 1/x as a function of T" is a straight 
line with intercept ^(q)°. 

Finally, combining the knowledge that for local moment systems, the magnetic response 
is governed by the ratio of the Zeeman energy gfisH to the thermal energy ksT, with eq.(^ 
yields another scaling law, namely 

ixiH, T)-i - e{Or/C)-'T" = f{H/T) (6) 

The law eq.(^ can be tested with extraordinary sensitivity in the long wavelength (g = 0) 
limit using bulk magnetometry. 

For systems in which the critical magnetic degrees of freedom derive from spin-density 
fluctuations of a Fermi sea, matters are much more complicated and do not satisfy eqs(|)-(|). 
Here we expect that non-local terms, derived from the band-structure, should have both E— 
and T— dependence. At low T, the first correction to the T = Pauli susceptibility is of 
order rather than T°, as suggested by the generalized Curie- Weiss Ansatz. Also, T and 
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iE are not interchangeable in the temperature-dependent susceptibihty. In particular, above 
two dimensions,0 the quantum critical fluctuations of a spin density wave are non-interacting 
at low energies and long wavelengths, in which case E/T scalingZl never occurs. Likewise, 
none of the less familiar forms@-(P) accounting for the q-, T-, x-, and if-dependence should 
obtain. Of special signiflcance is that an external magnetic fleld modifles the band structure, 
implying an influence on the non-local terms and an inability to scale the fleld-dependent 
effects simply in terms of the ratio g^sH /ksT^ 

We have thus reduced the problem of discriminating between the scenarios in flgs.[l|b) 
and c) to discovering whether the largely untested laws (^-(0) describe the T-, E-, 
and X- dependent magnetic response of CeCug-xAua^., and to verifying that with an order of 
magnitude improvement in energy range, E/T- scaling obtains at the critical composition 
and wavevector. To carry out this ambitious urogram, we have taken advantage of a low-T dc 
magnetometeril and the IRIS spectrometer^ at the ISIS spallation source. IRIS represents 
a substantial improvement over the instruments previously used to examine the magnetic 
fluctuations in CeCus.gAuo.i in that it (a) surveys the magnetic fluctuation spectrum over 
a very wide range of wavenumbers, (b) has an energy resolution 9/ieV = O.li^' almost an 
order of magnitude better than that of previous experiments, and (c) yields data in absolute 
units established using the simultaneously measured elastic incoherent scattering from the 
sample. 

As X is varied in CeCwe-x^^x, Bragg reflections due to antiferromagnetism appear at 
speciflc wavenumbers on a "butterfly" -shaped critical lineS, (black crosses in flg ||a) . At the 
quantum critical concentration x = Xc, the Bragg peaks are replaced by magnetic critical 
scattering - peaked on the same "butterfly"El. Fig. Qb) and c) show scans through the 
butterfly, obtained for the very small energy E = 35fieV. On warming from 50 mK to 
0.47^" ^ E/kB, the magnetic response at the butterfly is somewhat reduced, while on going 
to 1.5K = 4E/kB, it is greatly suppressed, showing that when T crosses E/ks, the magnetic 
response begins to change rapidly. Fig. ^ shows energy scans collected at 50mK at the 
various (nearly) flxed momentum transfers indicated in the inset and plotted following the 
new scaling paradigm eq(|^), using the exponent a ~ 0.75 obtained from the E/T scaling 
found in our previous critical wavenumber (q = (1.2,0,0)) experiment^ and verifled with 
much greater confldence below. The data collapse onto a single curve is clearly consistent 
with eq(|D, but there is sufficient scatter to mandate checks of the other untested laws (^ 
before we can claim to have verifled the local moment scenario giving rise to eq(|]). 

Fig.0 shows the modifled Curie- Weiss plots suggested by eq(|^) for CeCug-xAu^:, where 
the effective Weiss temperature is varied either by changing q in the neutron scattering 
experiment, or by changing x in the bulk (low fleld) magnetometry. Near the quantum 
critical point, the static susceptibility follows the modifled Curie- Weiss law over the entire 
two-decade temperature range of our measurements. As x moves away from Xc in either 
direction, x~^(T) becomes more curved as T — >^ 0, although there still is a quantum critical 
regime bounded by a cross-over temperature which increases with distance from Xc and above 
which the T" power law seems to hold (Fig. |^), and below which a tendency towards less 
singular (heavy) Fermi liquid behavior is seen. Fig.^a also demonstrates that at quantum 
criticality (i.e. x ~ 0.1) the same law is followed for general q, where is obtained 

via the Kramers-Kronig relation from the neutron data. Indeed, the parallelism of the lines 
in Fig.^ shows not only that a single exponent a is sufficient to account for all of the 
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magnetometry and neutron data, but also that the Curie constant C is - within the hmited 
statistical accuracy of the latter - the same for all wavenumbers. The overlap between 
the neutron data at q = (1.8,0,0) far from the butterfly but close to (2,0,0), equivalent to 
(0,0,0), and the q = bulk results is reassuring and shows that q = is a typical non-critical 
wavenumber, with a spectrum likely to resemble that measured at q = (1.8,0,0). A final 
result which can be extracted from Fig.^ is another estimate of the wavenumber- dependent 
Weiss temperature 6'(q), which is simply the x-axis intercept, taken to the power 1/a, of 
the lines going through the data for each q. The outcome is entirely consistent with that of 
the E/9{q) scaling analysis in fig.^). 

We finally turn to the energy- and T-dependent spectroscopies with the best signal to 
noise ratios, namely the neutron measurements of the scattering with q on the butterfly 
and the H- dependent static magnetization. Fig ^a shows the former, plotted versus E/T, 
collected at a wavenumber (0.8,0,0). This data set, which densely covers four decades of 
E/T, confirms that the scaling is optimal (see left inset) for a = 0.75 ± 0.05. Fig^ shows 
the collapse of the bulk magnetization data as suggested by the H/T scaling form (|^). The 
similarities in the scaling function and exponent {a = 3/4), obtained in this entirely inde- 
pendent assay, to those obtained via neutron measurements demonstrates that the magnetic 
response is driven by a singular local term. Two interesting details emerge. First, the ef- 
fective moment (for definition see Fig 4) = qhb = 1.5/iB is the same order as an atomic 
moment (0.6/iB)estimated from considering the g-space average of the neutron data inte- 
grated to E = ImeV. This rules out large, randomly occurring ferromagnetic clusters are 
a key feature of the quantum critical point.0 Second, the scaling function f{h) (where 
h = QfiBH/kBT) has the power-law asymptote h~°' at large h, which decays much more 
slowly than the exp{—h) asymptote derived from the Brillouin function. Thus, the active 
ingredient responsible for the critical behavior is an atomically local magnetic moment, with 
a size comparable with a single spin 5 = 1/2 with an intrinsically critical response to an 
external field. 

Our experiments reveal that the magnetic fluctuations in the heavy fermion metal 
CeCug-xAuj: are extraordinarily simple near the quantum critical point at x = 0.1. Not 
only do our new high-resolution neutron data satisfy E/T scaling at the critical wavenum- 
bers with much greater certainty than before, but they together with magnetization data 
satisfy two hitherto untested scaling laws and a generalized Curie- Weiss paradigm, all of 
which emerge from a local moment scenario for the quantum critical point. Detailed quan- 
titative analysis (see figure captions) reveals that an algebraic form. 
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C-' = Xo\E,T) + e{cir/C (7) 



accounts for all of our H- q- and T- dependent results. Eq(|^ is the most straightforward 
generalization of the Curie- Weiss description of a critical point and has only three inputs, 
an overall amplitude, an exponent a and the function 6{q). 

Clearly then, in answer to the question posed in the introduction, whether an itinerant 
(fig.|l|b) or local (fig.|l]c) picture of the magnetic quantum critical point applies to this heavy 
fermion system, it is the local scenario (fig.^c) that is appropriate. Beyond this, the most 
striking feature of the results is the appearance of a single unusual exponent a = 0.75 in 
the scaling properties at widely different wave-vectors. The E/T and H/T scaling with 
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the unusual exponent o; = 0.75 < 1 means that the magnetic properties of CeCue-ajAua; 
are more severely non- analytic than those of both Fermi liquids and conventional insulating 
magnets, where a = 1. 

We conclude that a Fermi hquid, CeCue, with one of the largest known effective car- 
rier masses is remarkably close to an instability marked by a new kind of critical behavior 
where the moments associated with the tightly bound f electrons are on the brink of separa- 
tion from the conduction electrons- how close is apparent from the approach of the inverse 
magnetic susceptibility to the anomalous T" asymptote as T rises above 2K( fig. 3b). The 
locality of the quantum fluctuations in the ordered phase suggests the development of a new 
kind of localized excitation, one derived from the ambiguity of whether the f-electrons are 
counted within the Fermi sea or not. These excitations arc neither heavy quasi-particles nor 
conventional spin excitations and might best be regarded as Fermi surface "shredders" : new 
degrees of freedom which are ultimately responsible for the non-Fermi liquid behavior and 
the unusual local scaling behavior at the quantum critical point. 
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Figure Captions 

Fig. 1. Phase diagram and scenarios for quantum criticality in CcCuq-xAux a) 
Phase diagram for CeCug-xAuj;, where increasing Au-concentration x drives the heavy 
fermion alloy into an antiferromagnetically ordered state. The quantum critical point 
at X = 0.1 is the subject of this paper. b)-c) Two competing models for the mag- 
netic quantum critical point in heavy fermions. W is the strength of the hybridization 
between the localized f-electrons and the surrounding conduction sea, which causes 
the antiferromagnetic order (with transition temperature T/v) to collapse at a critical 
value W = Wc- (b) The local moments are magnetically quenched - absorbed into the 
metallic bands- at a finite "Kondo" temperature, and do not play a role at the 
quantum critical point. Antiferromagnetism develops via a "spin-density- wave insta- 
bility" of the underlying Fermi surface, (c) Local moments exist down to an effective 
lattice Kondo temperature, T^, which vanishes at the quantum critical point. In this 
model, local moments become quenched precisely at the quantum critical point, and 
play an active role in the magnetic critical fluctuations. 

Fig. 2. Magnetic neutron scattering data, (a) Spectra at wavevectors q indicated in 
inset, as a function of energy transfer E in units of the q-dependent Weiss temperature 
6'(q). Intensities are correspondingly normalized to ^(q)~° (a = 0.75), as suggested 
by eq.(D. The line is the scaling function G{6) = C/(l + {-i6)") with 6 = E/{aT). 
The inset is a map of reciprocal space showing the critical line (black lines) close to 
the magnetic Bragg peaks in the ordered regime (black crosses) and the considered 
wave vector regions (labeled A-H). b,c, Wavevector (q) dependence of the magnetic 
fluctuations along the middle (b) and lower (c) trajectory in the inset of (a) for flxed 
E = 0.035mel^. x" = S{1 — e"^/'^^^) is shown, where the non-magnetic background, 
derived from scans with q parallel to the easy moment direction c, has been subtracted. 
Different colours represent different temperatures. Lines correspond to Eq.(7) by if = 
and a temperature independent 6'(q)" expanded in even powers of q. 

Fig. 3. Temperature dependence of the inverse static susceptibility, (a) x(q) (from 
direct magnetization (q = 0) and neutron data (q 7^ 0)), showing the same anomalous 
power-law dependence on temperature, over a wide range of different momentum values 
(see inset in flg.2a). b) Checking the same anomalous power-law in temperature of 
= 0) for different Au concentrations x. The temperature below which deviations 
occur increases with distance to the critical concentration x = 0.1. 

Fig. 4. Energy and field scaling Plots, (a) E/T scaling plots taken at various critical q 
vectors (labelled as A and B as in the inset of flg. (2a)), combining previous triple-axis 
data with the much more extensive time-of-flight data taken at IRIS. The line is the 
product of the Bose-Einstein factor (1 — e~^/^^^) multiplied by the scaling function 
g{y) = C/(l — iy)'^ with y = E/{aT) in eq.(3) with exponent a = 0.75. Right inset 
shows the wide range of E — T covered. Left inset shows the "scatter" of the scaling 
plot as a function of a, being minimal for a = 0.72 ± 0.05. (b) H/T-scaling of the 
local contribution to the uniform magnetization M{T,H). l/xo = {dM/dH)~^ — 
(4.1/i^/mev)^^. Solid line corresponds to scaling function f(h) = (1 + /i^)^°/^ with 
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h — guBH/ksT with a — 0.75 and an effective moment /i — g/iB — I-^/ib- Inset 
shows the H — T range where scahng apphes. 
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